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1. Introduction 

Event shapes measure geometrical properties of the energy flow in QCD final states. 
Their conceptual simplicity combined with a remarkable sensitivity to a range of 
features of QCD radiation, has led to their being among the most extensively studied 
QCD observables, both theoretically and experimentally. This is especially the case 
in e + e~ and DIS collisions, where seminal results include not just a plethora of 
measurements of the strong coupling |]], but also tests of the colour structure of 
QCD , detailed validation of Monte Carlo event generators || |, |^, H , and profound 
insights into the dynamics of hadronisation (for reviews see [0, |j ) . 



- 1 - 



Some of these studies consider mean values of event shapes, usually compared to 
fixed order perturbation theory, but the most information is to be gained from event- 
shape distributions. The majority of data corresponds to events whose energy flow 
closely resembles that of the lowest order (Born) event, a region usually associated 
with a small value (V) of the event shape. This is because large departures from the 
Born energy flow are caused by radiation of one or more hard gluons, each of which 
costs a power of a s . In contrast the dominant region, V <S 1, is characterised by the 
presence of multiple soft and collinear radiation and associated virtual corrections 
(Sudakov suppression), which leads to large perturbative terms (a$ In 2 V) n that have 
to be resummed to all orders. Fixed-order and resummed calculations, both at next- 
to-leading (NL) accuracy, have been vital to the full exploitation of the vast amount 
of data in e + e~ and DIS. 

At hadron colliders, event shapes (and a related class of observables, jet-res- 
olution threshold parameters) have so far received much less attention, there being 
for example at the Tevatron only few experimental measurements || |TD[. Recently 
however a number of tools for investigating them beyond leading order (or leading 
logarithms) have started to be developed JTT], 0, 0, [TJ, p], 0, [TJ, [1J. This 
opens up the possibility of using hadron-collider event shapes for a wide range of 
quantitative studies. 

Compared to e + e~ (and to some extent DIS) the hadron-collider experimental 
environment is more challenging, notably due to the presence of the underlying event 
(or equivalently beam fragmentation) and to the limited detector coverage in rapidity. 
Theoretically the situation is also more complex: while in e + e~ at lowest order there 
is radiation from just two outgoing jets, in hadron collider dijet events for example, 
at lowest order there is radiation from two incoming and two outgoing jets, with 
dependence on incoming parton distributions as well as a rich structure of interference 



between emissions from the different jets p0| , 21] , |22 , [23|, [24| , |2 



The more challenging environment and richer structure, though sources of tech- 
nical difficulties, provide many of the motivations for extending event-shape analyses 
to hadron colliders. Thus a study of hadron collider event shapes will not merely be a 
replication of what has been done in e + e~ and DIS (for example with measurements 
of a s ) but will allow novel investigations of both perturbative and non-perturbative 
aspects of QCD. 

Perturbatively, event-shape distributions are sensitive to the underlying jet- 
production channel: an event pp — > 2 jets with an underlying gg — ► gg partonic 
structure will typically have more radiation than qq — > qq, and this will be reflected 
in a distribution for that channel that is peaked at larger values of V. 

Dijet events will also allow tests of new structures in perturbative QCD — colour 
evolution matrices (also known as soft anomalous dimension matrices), initially de- 
veloped and studied at Stony Brook f2~I| , ^22] , |23| , |2^ , [25j — that arise for the first time 
in events with four or more hard partons. In processes with three jets (for example), 
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the combined colour charge of any pair of partons is uniquely determined, via colour 
conservation, by the colour charge of the third parton (triplet for a quark, octet for 
a gluon). This is important because the pair colour charge affects the pattern of 
large-angle soft radiation — the fact that it is unique means that one can ignore 
the evolution of colour associated with virtual corrections, allowing one to attribute 
a 'classical' probabilistic interpretation to the resummed virtual corrections. When 
there are four or more jets, colour conservation is insufficient to uniquely fix the 
combined colour charge of any given pair of partons, and one needs therefore to 
account for the 'quantum' evolution of colour in the QCD amplitudes, and this is 
accomplished with the Stony Brook colour evolution matrices. Some first investi- 
gations of their phenomenological consequences have been given in |24], |26| in the 
context of rapidity gaps, however event shapes should allow more extensive studies, 
as discussed later in the paper. 

On the non-perturbative front it is worth recalling that one of the most im- 
portant applications of e + e~ and DIS event-shape studies in recent years has been 
for extracting information about hadronisation, notably through 'power-correction' 
studies. It has been suggested ]27|, p8j , using renormalon @ inspired techniques, that 
leading hadronisation effects essentially act to shift event-shape distributions by an 
amount of order Aq C d/Q, with Q a hard scale of the process. The coefficient of 
this correction is given by the product of two factors: one of which (denoted Cy) is 
calculable in perturbation theory and observable dependent, the other one («o) being 
fundamentally non-perturbative, but universal for a range of observables in various 
processes. More sophisticated treatments, that attempt to go beyond leading effects, 
via 'shape functions' have also been developed P5| . There have been extensive tests 
of these ideas for 2-jet (qq) events, with remarkable successes both in e + e~ and DIS, 
as reviewed in ||. 

Many of the predictions from these renormalon-inspired techniques coincide with 
those from a simple Feynman-Field model of hadronisation p9| . To truly validate 
them requires that one investigate also processes with gluon jets and with non-trivial 
geometries of the jets (as in multi-jet events). Some studies in e + e~ and DIS are 
in progress, however multi-jet events there are rare (suppressed by powers of a s ). 
Hadron collider dijet events, in contrast, provide a natural environment for such 
tests, because of their inherent multi-jet structure and the large fraction of gluon 
jets. 

The beam jets do introduce certain complications in these studies, since their 
hadronisation (the part that is usually called 'underlying event') is not expected to be 
fully described by renormalon-inspired methods. Observables specifically intended 
for studying beam fragmentation have been proposed in |30| and examined in [|31 



Nevertheless, with the aid of a suitably complementary set of dijet event shapes, 
it should be possible to disentangle a number of features of the underlying event, 
in particular using approaches as discussed in fT2fl. In this respect it may also be 
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helpful to consider event shapes in the somewhat simpler case of Drell-Yan (with and 
without a jet). 

As in e + e~ and DIS, both fixed-order and resummed perturbative predictions 
at NL accuracy are essential in order to fully carry out these studies. Usually they 
are combined together via some form of 'matching' procedure. Fixed-order predic- 
tions for arbitrary infrared collinear safe observables can be obtained with the aid 
of fixed-order partonic Monte Carlos [JTSL O, [T5[. Resummed predictions, being 



more intricately linked with the details of the observable, had, till now, traditionally 
been obtained analytically, by hand. This made it tedious to consider more than a 
handful of observables, especially when dealing with the additional complications of 



multi-jet topologies 0, |32], [33|, [34], |35|. Recently however methods for automating 
resummations have been developed fl8| , [19], [36|] which enables a much wider range of 
studies. 

A current technical restriction with these automated resummations (as well as 
all fully NLL analytical multijet resummations carried out so far) is that they apply 
only to global event shapes p7| , i.e. observables that are sensitive to radiation in any 



direction, notably also close to the beam. 1 The only event shape distribution that 
has so far been measured at the Tevatron, a transverse thrust [R|, does not satisfy 
this property. 

One of the main purposes of this article is therefore to introduce and resum a 
set of global event shapes and jet-threshold resolution parameters. These include 
observables defined directly in terms of all particles in the event. Such definitions 
might seem particularly unsuited to practical uses, given that experimental detectors 
have limited angular reach close to the beam. We shall, however, discuss the impact 
on comparisons to theory of the limited experimentally available rapidity range, 
arguing that this is not an unsurmountable obstacle. 

We shall also introduce observables defined just in terms of particles in the 
central region of the detector, but for which one arranges an indirect sensitivity to 
the effects of all other remaining particles. Though these definitions eliminate the 
problems associated with limited detector reach, one should be aware that they are 
subject to certain extra theoretical (and potentially also experimental) problems. 

Throughout our discussions we shall refer to resummed results as obtained with 



the Computer Automated Expert Semi- Analytical Resummer (caesar) [fT9 |. We 
shall also consider observable-specific issues related to the sensitivity to the under- 
lying event. 



x It is to be noted also that the approximations (the angular ordering) contained in current 
Monte Carlo event generators such as Herwig || or Pythia Q are at their most accurate for global 
observables (Ariadne || in contrast has similar accuracy for global and non-global observables) — 
this means that global observables should be particularly suited also for the precise validation and 
tuning of event generators at hadron colliders. 
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2. General considerations 



In order to provide the background for the discussion of the various observables 
proposed later in this paper, it is helpful to first introduce some general concepts 
and notation. 

We shall consider here event shapes, which for any number N of final-state 
particles, are defined in terms of some function V(q 1 , . . . , q^) of the final-state four- 
momenta qi, . . . These event shapes will provide a continuous measure of the 
extent to which a general event's energy-momentum flow differs from that of a lowest 
order (Born) event. 

Typically one wishes to consider the event shape only for events that are suf- 
ficiently hard, requiring for example jets above some minimum transverse energy 
threshold E tjmin . We will denote this kind of hardness selection cut by a function 
H(qi, . . . , qiy), equal to 1 for events that pass the cuts and otherwise. One can then 
define the cross section for events that pass the cuts, 

<r w = W d$N^H( qi ,...,q N ), (2.1) 

where dajy/d^N is the differential cross section for producing N particles in some 
configuration While individual iV-particle cross sections are, of course, infrared- 
collinear (IRC) unsafe, the combined summed and integrated cross section is safe for 
any IRC safe function 7i. 

One also defines the partial integrated cross section T, n (v) for events that pass 
the cut and for which additionally the event shape observable is smaller than some 
value v, 

SwM = I Mn ^T Q ( v " V (*> • • • > *v)) H( qi , ...,q N ). (2.2) 

The differential normalised distribution for the event shape is then given by 

1 d^njv) 
a n dv 

2.1 Resummation 



(2.3) 



Resummations are relevant in the region of small v, where logarithmically enhanced 
contributions, (a s In 2 v ) n , are large at all orders, making fixed-order predictions un- 
reliable. The fact that the event shape is small implicitly means that the events 
resemble the Born event. This allows one to write a factorisation formula in which 
the resummed prediction for a given scattering channel 5 (for example qq — > qq, 
qq — > gg,...) T: n> s( v ) is expressed as the integral, over Born momentum config- 
urations £>, of the product of the differential Born cross section das/dB and the 
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resummed probability //3.5(f) that the event shape has a value smaller than v for 
that given Born configuration, 



Xh,s(v) = JdB^f BiS (v)H(p 3 ,p±), (2.4) 



where we use the convention that p 3 and P4 denote the momenta of the outgoing 
Born partons (while p\ , P2 are the incoming hard partons) . 

The resummed probability //3.5(f) can often be written in an exponential form 



fB,s( v ) = ex P [Lgi{a s L) + g 2 (a s L) + a s g 3 {a s L) H ] , L = In - , (2.5) 

where Lgi(a s L) resums leading- logarithmic (LL) terms, a"L n+1 , g 2 (a s L) resums 
next-to-leading logarithmic (NLL) terms, a"L n , and so on. Current state-of-the-art 
is NLL resummation. The exact form of the LL and NLL functions depends on the 
observable under consideration. 

One can characterise a given observable by its functional dependence on the 
momentum of a single soft emission, collinear to one of the hard ('Born') partons in 
the event. For all known event shapes, this can be written as 

V({p},k) = dA-^j e-^%(0), (2.6) 

where {p} denotes the Born momenta (including recoil effects) and k is the soft 
collinear emission; fcj and 77 ^ denote respectively its transverse momentum and 
rapidity, as measured with respect to the Born parton ('leg') labelled l\ is the 
azimuthal angle of the emission with respect to a suitably defined event plane (when 
relevant); and Q is the hard scale of the problem. For hadron-collider dijet events 
we will use the convention I = 1,2 for the incoming hard partons and I = 3,4 for 
the outgoing hard partons. 

The values of the coefficients ag, bg and the form of the function gg(4>) in eq. (|2.6j ) 
are among the main characteristics of the observable that enter in the resummation 
formula eq. ( |2.5|) . For example the leading logarithmic function Lgi(a s L) depends 
only on the ag and bg values, as can be illustrated from its expansion 



L 9l (a s L) = - £ . — + O (a s 2 L 3 ) , (2.7) 

^— ' ag(a e + bg) it 



where Cg is the colour charge (Cp or Ca) of hard parton i. For the observables we 
consider here — of the continuously global |37|, |I| variety — the ag are all equal, 
di — ct>2 — ■ ■ ■ = a (and the dg are all non-zero). Because of the appearance of the 
hard-parton colour factors in eq. ( |2.7| ), the LL terms depend also on the underlying 
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hard-scattering channel 5 (e.g. qq — > qq as opposed to gg —>■ gg), though not on the 
particular momentum configuration B of the hard partons. 

The NLL terms are somewhat more complex. One can often separate them into 
two pieces, 

g 2 (a s L) = g 2s (a s L) + \nF(R'(a s L)) . (2.8) 

The first term g 2s (a s L) depends just on the 'single-emission' parameters of eq. ( |2.6|) 
(including de and ge{4>)), and on the Born momenta B and the channel 5. It contains, 
among other contributions, the dependence on the Stony Brook colour evolution 
matrices that were mentioned in the introduction. 

The second term in eq. (|2.8|) , !F(R'), is the only part of the NLL resummed 
formula that is affected by the observable's dependence on multiple emissions. It is 
a function of 

R'(a s L) = -d L L gi (a s L). (2.9) 

Schematically, for R 1 <C 1 (L not too large, v moderately small) there is one soft and 
collinear emission in the event that is harder than all other emissions. This hardest 
emission gives the dominant contribution to the value of the observable 2 and only the 
single-emission properties of the observable are relevant, so T(R') — 1. For R' ~ 1 
[L ~ l/a s , ?; < 1) there are typically R' + 1 similarly hard gluons that dominate the 
structure of the event. If multiple emissions tend to increase the value of the event 
shape, then for a given value of the event shape, this must be compensated for by an 
extra suppression of emissions, i.e. F(R') < 1. Conversely if the effects of multiple 
emissions tend to cancel out, this is compensated by reducing the overall amount 
of suppression, T(R') > 1. The function jF(i?') often depends on the underlying 
scattering channel S, but not on the hard momentum configuration. 

Throughout this paper, much of the discussion of the properties of observables 
will be framed in terms of the values of the ae, be, de and the forms of ge{4>) and 
F{R'). 

2.2 Experimental issues 



Our understanding of the current Tevatron detectors pl| , [f2|, [43], |44j] is that typically 
they have good measurement capabilities in a central detector region, up to about 
2.5-3.0 units of rapidity i], as well as a reasonable degree of measurement capability 
in the semi-forward region up to a rapidity of about 3.5, with much more limited 
capabilities up to a rapidity of about 4.5-5. The LHC detectors EH HB| should 
provide good coverage up to about 5 units of rapidity. 

A limit on the maximum reach in rapidity can be a source of problems because 
the most natural way of defining global observables is directly in terms of all particle 
momenta in the event. Introducing a cut on the maximum rapidity at which one 



2 Strictly speaking this is only true for recursively infrared-collinear safe observables [|l8| |l9| 
most event shapes are in this class. 
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measures particles potentially creates a mismatch between the 'ideal' theoretical 
definition and the measurement. 



It was pointed out however, in [|12j| , that, for sufficiently large values of the 
maximum accessible rapidity, r? max , the excluded kinematic region gives at most a 
small contribution to the observable, which for the discussion of a general observable, 
eq. ( |2.6| ), corresponds to an amount of order e _ ( a+6min ^ max where b min = min{&i,&2}- 
Accordingly, as long as one is in a region where the observable's value is larger than 
this, equivalently 

L < (a + b min )i] mSLX , (2.10) 

the resummed prediction should remain valid. One of the features of the resummed 
results that we shall therefore examine closely in the subsequent sections, is the 
typical values taken by the observable. As long as the cross section outside the 
region eq. ( 2.10|) is small, the effect of the limited experimental rapidity reach should 



be negligible. As a reference value for the maximum attainable rapidity, to be used 
in eq. (|2.10| ) when comparing with actual resummed results, we shall take ^ max = 3.5. 



A point worth mentioning is that measurement of the contribution to an event 
shape from the forward region may be somewhat simpler than detailed jet studies 
in that region. This is because many event shapes' sensitivity to the forward region 
is simply in terms of the total transverse energy deposited there, so that it is not 
necessary to resolve the detailed r/-0 structure of the energy deposition. 

An alternative to defining global observables in terms of all momenta in the event, 
to be discussed in section ||, is to use particles only in some restricted central region 
C while introducing an indirect sensitivity to momenta outside that region through 
the addition of the total vector sum of transverse energy in C, or equivalently, the 
missing transverse energy. In such cases, the critical experimental issue will no 
longer be the rapidity reach of the detector, but rather the accuracy with which the 
missing transverse energy can be determined. If the missing transverse energy can 
be determined to within an error 5E±, then the accessible region for studying the 
event shape will be roughly L < a\n(E±j et /5E^), where E±j et is the jet transverse 
energy. 



2.3 Event selection cuts 

The selection cuts that we propose for event-shape studies are as follows. Using some 
IRC safe jet-algorithm, one should first select events where the two jets with highest 
transverse energy (E± t x > E± <2 > • • •) are both central, | ??jets 1,2 1 < Vc, where the limit 
of the central region i] c is taken to be of the order of about 0.5 to 1. 

Restricting the jets to the central region is advantageous theoretically, because it 
eliminates large subleading corrections that are likely to be associated with forward 
jets. It should also be advantageous experimentally, insofar as the measurement of 
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the event shape needs particularly fine resolution in rapidity and azimuth in the 
region containing the jets. 

A cut should also be placed on the hardness of the jets. To avoid problems of 
perturbative convergence of an that arise with the use of symmetric E± cuts on the 
jets j|7|, fffi| , one can for example place a cut on the transverse energy E± ) i of the 
hardest jet, E ±Mn < E ±jl < E ±>max . 

At the (NLL) accuracy that we discuss in this paper, the details of the jet algo- 
rithm used to identify the two hardest jets, be it of cone [[49] or cluster |50], type, 
do not modify the predictions. This is a consequence of the continuous globalness 
of the observables and the fact that we consider just the soft and collinear limit for 
the observable. At higher logarithmic accuracy, there are reasons to believe that 
the partial integrated cross section £^(t>) remains independent of the jet algorithm. 
However since the total selected cross section, an, does depend on the jet algorithm 
at NLO, so does the ratio eq. ( |2.3| ) at NNLL. In the region of finite v, ^n( v ) itself 
also depends on the jet algorithm. 

In the plots shown throughout this paper, we take r\ c = 0.7, E± jto i n = 50 GeV 
and place no limit of -E_i_ jma x- We assume a Tevatron run II regime, pp-collisions 
at the centre-of-mass energy y/s = 1.96TeV and use the CTEQ6M parton density 
set p^| , corresponding to a s (Mz) = 0.118. 

We set the factorisation and renormalisation scales to be the sum of the trans- 
verse energies of the two most energetic jets. Future work |53[ will include fixed-order 
matching and will explore a range of alternative scales and perform a systematic 
study of the dependence on these scales. 



3. Directly global observables 

We will first consider observables that are defined in terms of all emissions in the 
event. This direct sensitivity to all emission momenta is the origin of the name 'di- 
rectly' global. There will be two groups of observables — those that can naturally be 
defined in terms of all emissions in the event, and those which would more naturally 
be defined in terms of emissions in some central region, to which we add an extra 
term that is sensitive to forward emissions. In this section we will consider the first 
class of observables, while the second kind will be discussed in section |. 

3.1 Transverse thrust (a detailed example) 

We give here a detailed discussion of the 'transverse thrust' event shape. Since the 
issues that arise are quite similar for nearly all event shapes, the discussion for the 
other event shapes that we define afterwards will be somewhat briefer, restricted to 
highlighting important differences compared to the transverse thrust. 
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The most obvious extension of the various e + e event shapes is to define a thrust 



>g in the transverse plane, 



max ■ 



(3.1) 



where the sum runs over all particles ^ in the final state, q±i represents the two 
momentum components transverse to the beam 3 and Ht is the transverse vector that 
maximises the projection. 

Variants of the transverse thrust based just on particles in a restricted central 
region have been measured in [llll] and calculated at fixed order in Those 



] and calculated at fixed order in 
variants are discontinuously global (a 12 7^ ^3,4), and so beyond the scope of the 
current automated resummation technology. 
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Table 1: Leg properties for tj_ j9 . Here G ~ 0.915965594 is Catalan's constant. Starred 
entries (here and in subsequent tables) indicate quantities that have been determined only 
numerically by CAESAR, and for which we have manually provided the full analytical in- 
formation. 

The transverse thrust as defined in 
eq. ([34]) has already been used as brief 
example of the capabilities of CAESAR 
in |Tg|] , but here we discuss it more ex- 



tensively. In the 2 + 2 jet limit it tends 
to T±„ = 1. Accordingly the observable 



CD 



1—T] g . Its de- 




Dependence 



Figure 1 

the azimuthal angle 
collinear to each leg 



of T±„ on 



that one resums is r± tS 
pendence on soft and collinear emissions 
is given by eq. ( f2.6|) with the parameters 
shown in table [3]. The fact that bt = 
for the incoming legs £=1,2 means that 
along the beam direction the observable 
is sensitive to the emitted transverse mo- 
mentum uniformly at all rapidities. In 
contrast for the outgoing legs £ = 3,4, 

we have be = 1, meaning that the effect of emissions close to the legs is suppressed. 
From eq. (|2.7|) , one sees that this translates to smaller LL contributions from the 
outgoing jets than for the incoming jets. 

3 We distinguish between which always refers to a transverse momentum with respect to the 



of a single emission 



beam direction, and H', used in eq. (2.6), which is a transverse momentum with respect to a given 
leg. 
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The observable's dependence on the azimuthal angle of a single emission, ge(<fi), 
is shown in fig. [I]. For both incoming and outgoing legs, the observable is most 
sensitive to radiation perpendicular to the event plane, = | , The analytical 
form for gi{4>) has been established by CAESAR only for the outgoing legs, it being 
given numerically for the incoming legs (though the analytical form is straightforward 
to derive by hand, and reads #1,2(0) = 1 — I cos0|). 

In order to get an understanding of the coefficients di, one needs to know the 
value of the hard scale Q in (|2.6| ). For the present analysis (and for all other analyses 
in this paper) we have chosen to set this scale to the partonic centre of mass energy of 
the hard collision\/5 and have considered a hard configuration in which the scatter- 
ing angle 9 between the outgoing jets and the beam direction (in the hard-scattering 
centre of mass) satisfies cos# = 0.2. CAESAR would provide in table [3] the numer- 
ical value of the coefficient di corresponding to this reference hard configuration. 
However, for all the observables considered in this paper, we report the functional 
dependence of dt on the angle 8, as well as the explicit form of gi{4>), since they 
can be derived from simple analytical considerations (or by examining the numerical 
results provided by caesar). To distinguish then the genuine output of CAESAR 
from the quantities obtained by hand, the latter are highlighted with an asterisk. 

Table [3] also includes the result for (ln^(0)) = J ^Llnge^). This is of interest 
insofar as it is actually the combination lnci^ = lnci^ + (In ge(4>)) that appears in the 
resummation formulae, rather than di or ge(<p) separately. 

In order to complete the information needed for the resummation, caesar evalu- 
ates also the function jF(i?') that appears in eq. (|2lf). The transverse thrust (strictly, 
T ±,g) has the property of additivity, meaning that in the presence of many emissions, 
the value of the observable is simply the sum of the values that the observable would 
take for each emission individually. For observables with this property, T{R') is 
known analytically, 

p -7eH' 

This completes the information needed so as to make resummed predictions for the 
transverse thrust. 

The resulting differential cross sections for different partonic scattering channels 
are shown in figure ^|, where q denotes a generic fermion (quark or antiquark) and 
qg — > qg indicates any process with incoming fermion and gluon. They have been 
obtained by integrating over all hard events that satisfy the cuts discussed in sec- 
tion |2l| The dg are redetermined for each hard configuration, since they depend on 
the kinematics of the event. While in table [3] the arbitrary reference hard scale was 
taken to be v^§, hereon, for the purposes of calculating differential cross sections we 
choose Q, as well as the renormalisation and factorisation scales, to be the sum of 
the two hardest jet transverse energies, which is closer to the virtuality of the ex- 
changed parton in the hard scattering. Since the distributions are purely resummed, 
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Figure 2: The differential cross section 
<iS^/dlnr_|_ 9 separated into the different 
hard-scattering channels. 
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without the inclusion of the exact fixed-order contribution (matching, left to future 
work they are reliable only at small values of t± j9 . For t_i_ i9 ~ 1, the differential 
cross sections are negative, this being a reflection of a breakdown of the resummation 
approach in that region. 

At small T± t9 the differential cross sections of fig. |2| have somewhat different 
shapes for the various hard-scattering channels (S). This is a consequence of the 
different colour factors multiplying the LL terms in eq. fl2.71) , the LL contribution 
being larger for channels with more hard gluons — such channels tend to radiate 
more, so the distribution is peaked at larger values of T± >g . 

Combining the different hard-scattering channels, one obtains the full resummed 
distribution, fig. ||. In section |2]2] we discussed the consequences of limited experi- 



mental rapidity coverage for comparisons between (fully global) theory and (partially 
global) data. The main result was that the resummation remains valid, in the range 
limited by eq. (|2.10|) , which here translates to L < r] max . A subtlety arises in the 



translation between t_i_ iS and the value of L: as discussed in appendix |A], in eq. (fX 
we resum not lnl/fi^ but actually In l/(Xr± >g ), where X is chosen so as to cancel 



the average value of dege((J)). For the hard events satisfying the cuts in section |2.3| this 



corresponds to X ~ 5. This is a convention, adopted from e + e~ [53| and DIS ffl], 
intended to minimise spurious subleading logarithms associated with the potentially 
arbitrary normalisation of one's observable. 4 Since for the discussion of section \2.2[ 



4 Readers familiar also with |18| may have noticed that the cross sections shown in fig. 1 (there) 
and fig. [2] (here) differ considerably in normalisation, even though they apply to the same observable. 
This is because in [|l8) we used different cuts (r] c = 1.0 rather than 0.7), but also because there we 
used X = 1 and Q = v§, both choices being associated, for this observable, with large subleading 
logarithms. 
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it is L = In l/(Xr± >g ) that determines the maximum value of rapidity that is relevant 
in the resummation, eq. ( [2.10|) should be rewritten 



\nv>-(a + b £ )(r] max + {\nX)). (3.3) 



Assuming r] max ~ 3.5, one obtains the limit In r^ iS > —5, suggesting that the distribu- 
tion is nearly entirely in a region where the limited experimental rapidity reach should 
be unimportant. If however one goes to higher energies, such as E± t \ > 200 GeV as 
shown in fig. ^|, channels with smaller colour factors start to dominate (e.g. qq — > qq), 
because one samples the parton distributions at higher x; this together with the 
smaller value of a s leads to the r± i9 distribution being dominated by lower values of 
the observable. 

A final point regarding t±„ concerns its sensitivity to the underlying event. To 



parametrise this sensitivity it is useful, as was done in f!2 |, to take a simple model 
[p0|1 in which particles from the underlying event have a spectrum dn^ u ' e ^ / dk±dr] 
that is independent of r] and of the underlying hard scattering process. Since r_i_ i9 is 
additive, and sensitive along the incoming legs to (see tab. [1]), it receives a 

mean contribution from the underlying event 



{S r± ' 9> = U*_ ' + ( 7 J dk± ^ dHSklE^TEZ 2 (3 ' 4a) 

= < *Py > (2 w + 0(l)), (3.4b) 



where (k^' e ^) is the mean transverse momentum per unit rapidity coming from the 
underlying event. The unspecified contribution of O (1) reflects the fact that for this 
qualitative discussion, we have not taken into account the details of the observable's 
sensitivity to emissions close to the outgoing jets. Though the result is for the mean 
effect of the underlying event, expectations based on other studies of non-perturbative 
effects |j5] suggest that the distribution of T± id will simply be shifted by (<^ u ' e 'Vj_ iff } . 
Because of the proportionality to 77 max , the effect will be large, making tj_ )S a good 
observable for testing models of the underlying event. For example the effects of the 
underlying event could depend on the underlying hard scattering channel, 5 and it 
might then be possible, from the event-shape distribution, to establish whether or 
not this is the case. 



3 One of us (GPS) wishes to thank B. R. Webber for discussions on this point. 
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3.2 Thrust minor 

Given the transverse thrust axis Ht, one can define a directly global thrust minor, 6 



T = 

-*■ m.g — 



i Qxi 



(3.5) 



where the x direction is defined as that perpendicular to the beam and to the global 
transverse thrust axis, which together define the event plane. The thrust minor can 
be viewed as a measure of the out-of-plane momentum, and is quite similar to the 
observable resummed analytically in [12] for Drell-Yan plus jet production. 

One can also normalise the thrust minor to E± tl + E± j2 - This will not modify 
the resummation, but it will affect the distribution at values of T ms of order 1. 
(Such a change cannot be made for the Tj_ jff , because it would destroy the positive 
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Table 2 

(ln^(0)), not shown, is 



Leg properties for T m ,g\ 



In 2 for all legs. 



definiteness of r± t . 

The main difference between the thrust 
minor and the transverse thrust concerns 
the properties of the outgoing legs, which 
for T m 9 have b £ = (table 0). Referring to 
eq. ( |2.7| ), one sees that this implies larger 
double logarithms than for t± j9 . Addition- 
ally, the dependence on multiple emissions 
is more complex, T mj9 not having the sim- 
ple property of additivity, so it is necessary 

to calculate jF(i?') numerically. The result, shown in fig. |4], is similar (though not 
identical) for all hard-scattering channels, and uniformly below 1, which is indicative 
of the fact that multiple emissions increase the value of the observable relative to 
a single emission. We also note that the function JF does not depend on the Born 
momenta £>, as can be verified numerically. 

The larger double logarithms (63^ = 0) and the larger coefficient di (together 
with the larger average value of ge(4>)) all contribute to T mg having a distribution, 
fig. ||, that is dominated by considerably larger values of the observable than was the 
case for r± t9 . As a result the peak of the distribution is actually at only moderately 
small values of T mjff , where fixed-order corrections (matching) may have a significant 
effect on the distribution. 

The rescaling factor, X, in the resummation here (of lnl/(XT m ^ g )) is X = 1. 
As a result the limited detector reach translates to a limit of applicability of the 
resummation lnT mj9 > —r] max . The only channel that extends at all below our 
reference value of 77 max — 3.5 is qq — > qq, which however, at these energies does not 
contribute significantly to the total distribution. 

Finally we observe that the linearity of the definition of T mg implies that it 
should have a sensitivity to the underlying event that is qualitatively quite similar 
to that of 7j g . eqs. (|3~4j). 



6 In the literature ||, JOj, such an observable, with an additional cut on the rapidities of measured 
particles, cf. eq. (4.5), has been referred to as a broadening. Our choice of nomenclature is in analogy 
with the thrust minor as defined in e + e _ [dm. 
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3.3 Three-jet resolution threshold 



Clustering jet algorithms, such as the fcf-algorithm |50|, |5T], |>7]], involve a successive 
set of pairwise particle recombinations, so as to group particles into jets. These 
algorithms often involve a jet-resolution parameter, which determines the point at 
which to stop the pairwise recombination. The smaller the parameter, the larger the 
number of separate jets that gets resolved. 

One of the most widely studied observables in e + e~ is the jet resolution threshold, 
t/23, above which the event is classified as having two jets, and below which it has three 
(or more) jets. Though an analogous threshold has been defined for going from 2 + 2 
to 2 + 3 jet events with the longitudinally invariant /c r algorithm in hadron-hadron 
collisions [ 50f1 , it has received somewhat less attention. Resummed calculations [11| 



(to lower accuracy than is given here) and measurements [|58] have been carried out 
only for the related problem of subjets within a single jet. 

Jet threshold resolutions and jet rates are closely related since the integrated 
cross section £«(?/) for y 2 3 is equivalent to the cross section for having a two-jet 
event, given a jet resolution of y. 

We will consider here just one variant of the fc r algorithm (the exclusive form of 
that adopted for run II of the Tevatron Reference implementations of this and 

other schemes are available in [55 . 



1. One defines, for all n final-state (pseudo) particles still in the event, 

d kB = ql k , (3.6) 
and for each pair of final state particles 

d M = mm{q 2 ±k , q\i\ ((% - Vif + {<f>k ~ 0/) 2 ) • (3-7) 



- 15 - 



2. One determines the minimum over k and I of the du and the c4b and calls 
it dS n >. If the smallest value is diB then particle qi is included in the beam 
and eliminated from the final state particles. If the smallest value is dij then 
particles q^ and qj are recombined into a pseudoparticle (jet). A number of 
recombination procedures exist. We adopt the E-scheme, in which the particle 
four-momenta are simply added together, 

Qij = Qi + Qj ■ (3.8) 



3. The procedure is repeated until only 3 pseudoparticles are left in the final state. 
The observable we resum is then 



2/23 = imax{rf (n) } , 



(3.9) 



where E± is defined by further clustering the event until only two jets remain 
and taking E± as the sum of the two jet transverse energies, 



E± = E±i + Ej 



(3.10) 



The reason for considering max n > 3 {d^} in eq. ( |3.9| ), instead of simply d^ 3 \ is 
that the recombination procedure is not necessarily monotonic in the d^ [|50|1 , 
so that n > m does not automatically imply <i <n - ) < d^ m \ 
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Table 3: Leg properties for 1/23. 



From table |3| one sees that, for all legs, 
the effect of single emissions scales as the 
squared transverse momentum, without any 
rapidity or azimuthal dependence — in this 
respect the observable is similar to 7/23 i n 
e + e _ . 

An interesting theoretical feature of the 
hadronic 2/23 concerns J 7 , fig. |6|. In contrast 

to the transverse thrust and thrust minor cases, T{R') has a strong dependence on 
the scattering channel. It seems that along the outgoing legs, multiple subjets can 
combine together before being combined with the main hard jet — this means that 
multiple emissions along an outgoing leg lead to a larger y 2 3 than a single emission, 
causing T to be smaller than 1, as in e + e _ [36]. For incoming legs instead, subjets 
tend not to recombine together, so that the value of 7/23 is determined exclusively 
by the hardest subjet. This type of situation leads to T{K) = 1. The actual value 
for T then depends on the relative importance (colour factors) of the incoming and 
outgoing jets. 

Figure shows the distribution of 2/23 for different hard subprocesses. It extends 
to considerably smaller values than, say, T mt9 , the reason being that, from the point 
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of view of LL contributions 2/23 ~ ^m,g- This suggests that the distribution of \ In 2/23 
should be comparable to that of lnT mg . One sees that to some extent this is the 
case, though the 7/23 distribution is somewhat wider than would be expected based 
on this argument. This is a consequence of the larger suppression that T m 9 receives 
from its J-'(R') (note that at a given value of y 23 ~ T^ g , the definition of R' is such 
that R' ma = ii^ m J. 

The fact that 2/23 ~ is important also from the point of view of the impact 
of the limited experimental rapidity reach, since the range in which the cut can be 
ignored theoretically is lny 23 > — 2^ ma x- One sees that for ^ max = 3.5, this includes 
most of the distribution, except in the qq — > qq channel, which extends a little 
beyond. 

While perturbatively 1/23 and T^ g are similar, they should have very different 
sensitivities to non-perturbative physics. It is already known that in e + e~, fct-type 
jet algorithms are relatively insensitive to non-perturbative effects. The same should 
be true also of the hadronic 1/23 — for example, the fact that along incoming legs it 
is the hardest sub jet that dominates, means that unlike T mjS , y 2 3 does not 'collect' 
non-perturbative contributions from the whole range of 1 77 1 < r/ max , but rather just 
from a limited region around the hardest subjet. 



4. Observables with exponentially suppressed forward terms 

While the observables given above could be easily defined in a 'global' manner, other 
jet observables, such as jet masses or broadenings, are more naturally defined just 
in terms of a central region, or of the two leading jets. As a result they are usu- 
ally not global, which prevents them from being resummed with current automated 
approaches such as CAESAR. 
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One solution to this problem is to explicitly add to them a term sensitive to 
emissions along the beam direction, so as to render them global. If we add just the 
q±i, we will obtain observables whose incoming legs have rather similar properties 
to the transverse thrust or the thrust minor, in particular with regards to their 
sensitivity to the underlying event and to relevance of the experimental ?7 max cutoff. 

Instead we add a contribution, ~ q±iC~^, with an exponential suppression in 
the forward direction. 7 This will usually lead to observables having bi = a on the 
incoming legs, which doubles the range of L, eq. ( p.lOj ), in which the resummed 
predictions are insensitive to r] max . Such observables will usually also have a reduced 
sensitivity to the underlying event, since the integral eq. ( |3.4a|) will converge rapidly 
in r] giving a result of order (k^' e ^)/ (E± t i + E± :2 ), without any ?7 max enhancement. 

Before defining such observables, let us first discuss the central region, C, in 
which the 'main' event shape (for example the sum of jet masses) will be measured. 
Various options are possible, for example taking all particles in the two hardest jets; 
or taking all particles having \rji\ < t] c + 5t], where r] c specifies the region in which the 
two hardest jets should lie (section pT3| ) , while 5r] extends this region, so as to ensure 
that all selected jets are well contained within C. For implementation in CAESAR, we 
have used this second definition, with 5r] = 0.4, though it should be noted that since 
our final variables will be global, the NLL resummed results are actually independent 
of the precise definition of C. 8 

Having established C, we then introduce the mean transverse-energy weighted 
rapidity, rjc, of the central region, 

Vc = y^—y^ViQi-i, Q-L.C = J] QUi , (4- 1 ) 

and define the exponentially suppressed forward term, 

£<J = q±i e~ lril ~ vcl , (4.2) 

such that it is invariant with respect to longitudinal boosts (modulo the dependence 
of C itself on boosts). We are now in a position to define some observables. 

4.1 Transverse thrust, minor and jet resolution 



Transverse thrust. Let us first define a thrust similar to that of section |3J], but 
only in terms of the particles in C, 



Tj. iC = max , t x ,c = 1 - T± tC ■ (4.3) 

n T,C V±,C 



7 In the spirit of [[30| |6l[] one could also investigate continuous classes of observables, in which 
the forward term goes as gj_ie _c ' ,K ' , c being a parameter that allows one to span a whole class of 
observables. 

8 While at higher resummed orders, and in the fixed-order perturbative and the non-perturbative 
contributions there will be dependence on the definition of C. 
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We then add the contribution £q with 
the dependence on the forward emission, 

T±,£ = T ± ,c + £c, (4.4) 

to obtain a global observable, 9 the resulting 
leg properties being shown in table |]. One 
sees that for the outgoing legs, they coincide Table 4 . Leg properties for T± g 
with those of table [l] for Tj_ )ff , while for the incoming legs — 1, as expected. Like 
T|_ iff , the observable is additive, so T is known analytically. 

We could at this point continue the discussion for T± t g along lines similar to those 
for the observables given above. Since, however, one can imagine introducing quite 
a few further observables, we prefer from now on to concentrate on just a subset of 
them, so as to illustrate interesting new features. Accordingly we refer the reader 
to detailed web pages [62] for further information about T± t £, as well as analogous 
extensions of the thrust minor, 

1 



- m,C 



,c 



T 



m,£ 



-m,C 



C ! 



(4.5) 



iec 



(because of kinematic recoil, it turns out that T m £ is global, and T m £ has &i )2 = and 
identical resummation properties to T m> c); and the three-jet resolution threshold, 
with |/23,e defined by the algorithm of section |0| applied only to the final state 
particles in C (and the beam) and 



V23,£ = V23,C + £j 



(4.6) 



Note that it is necessary to add £j here (rather than Sq), so as to ensure the contin- 
uous globalness of the observable, a\ = a 2 = a 3 = a 4 . 

Let us now concentrate in more detail on some observables that are more natu- 
rally defined without explicit reference to a central region. 

4.2 Jet masses 

Having determined a (central) transverse thrust axis nr,c as above, one can separate 
the central region C into an up part Cy consisting of all particles in C with p±-nT,c > 
and a down part Co, particles in C with p±_ ■ nxfi < 0. If C is taken to be made up 
of all particles in the two hardest jets, one can also define Cu and Co as consisting 
of the two jets separately. Such an alternative definition should not change any of 
the NLL resummed predictions for the global variables, and may actually help to 
minimise subleading (NNLL) corr ect ions . 10 

9 It actually turns out that even without the addition of any direct dependence on particles not 
in C, ti.c has an indirect sensitivity to them via effects of recoil of the hard jets. This leads to rx,c 
being discontinuously global (di,2 = 2, 03 4 = 1) rather than non-global, which is still beyond the 

SCOpe of CAESAR. 

10 In contrast, for the central, non-global observables (those with a C suffix), the exact choice of 
C will affect the NLL terms. 
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One then defines, in analogy with e + e [S3], the normalised squared invariant 
masses of the two regions 11 

2 
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(4.7) 
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(4.9) 
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Table 5: Leg properties for ps,e (and 
PH,e)- 



from which one can obtain a (non-global) central sum of masses and heavy-mass, 

ps,c = pu,c + Pd,c , Ph,c = max{p t / iC , p Dfi } , (4.8) 

together with versions that include the addition of the exponentially-suppressed for- 
ward term, 

PS,£ = PSfi + £■£ > PH,£ 

The single-emission leg properties for p$,£ 
and ph,e are identical (even beyond the soft- 
collinear limit), and quite similar to those 
of T ±,e except for the lack of azimuthal de- 
pendence on the outgoing legs. The sum of 
squared masses, ps t e, is an additive observ- 
able, so J-'(R') is given by eq. (|3.2|), while for 
the heavy-mass it needs to be computed nu- 
merically (and decreases for increasing R'). 

The resummed distribution for p H£ is 
shown, channel-by-channel in fig. Com- 
pared, say, to an observable like T mt9 , the 
distributions extend to relatively low val- 
ues, a consequence of the smaller LL terms 
(since h = a). However, as mentioned at 
the beginning of this section, the fact that 
b\ t 2 — a — 1 means that the limit of appli- 
cability of the resummation, in the presence 
of an experimental rapidity cut, is given by 
m PH,£ > — 2r? max . Thus all channels for 
the distribution are comfortably contained 
within this region. Figure 8: The distribution of pn,£ in the 

different hard-scattering channels. 

4.3 Jet Broadenings 

With the same division into up and down regions as for the jet masses, one can define 
jet broadenings. To do so in a boost-invariant manner, one first introduces rapidities 
and azimuthal angles of axes for the up and down regions, 
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(4.10) 



11 For certain non-perturbative studies, it can be advantageous [|64| to use a so-called 'E'-scheme' 
definition in which the 3-momenta ql arc rcscalcd ql — > (Ei/\qi\)qi, equivalent to a massless approx- 
imation for all particles. 
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and defines broadenings for the two regions, 
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(4.11) 



from which one can obtain central total and wide-jet broadenings, 

B T ,c = B Ufi + B Dfi , B Wfi = max{B Ufi , B Dfi } . (4.12) 

Adding the exponentially-suppressed forward terms gives global observables, 
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(4.13) 
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The single-emission properties of these two 
observables are identical and are shown in ta- 
ble |6|. In both cases J-'(R') needs to be cal- 
culated numerically, and it decreases with in- 
creasing R' (more strongly for -Br.f)- 

The distribution for Bw,e, separated into 
channels, is shown in fig. |9|. As expected from Table 6: Leg properties for Bt,s (and 
the different double logarithmic structure for Bw,s)- 
the outgoing legs, the B W £ distributed is cen- 
tred at larger values of the observable than pu,s- This can be seen clearly also from 
fig. ITU] which shows both distributions, summed over channels for two different E \ )Tn in 
cuts. Since a limited experimental rapidity translates into an allowed range of the 
observable's value \nv > — 27/ max , we see that the distributions for both p H£ and 
By\/£ are within the theoretically reliable region. The same statement applies to 
most of the other observables defined in this section. 



5. Indirectly global observables 

We call 'indirectly global' those observables that explicitly measure only a subset of 
the particles in an event, but are nevertheless indirectly sensitive to the remaining 
emissions, typically through recoil. Previously studied examples include certain DIS 
Breit-frame event shapes, such as the current-hemisphere broadening with respect 
to the photon axis or one of the variants of the thrust, t z q |pHfl . 

Indirectly global observables have the advantage that they can be defined purely 
in terms of particles in the central region C, while still remaining global. This elimi- 
nates potential problems associated with the limited experimental reach in rapidity. 

To construct an indirectly global observable, one takes one of the 'central', non- 
global observables of the previous section, and then adds to it some second 'recoil' 
quantity, also defined only in terms of the momenta of C, but that is sensitive to 
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the emissions outside C. An example of such a 'recoil' quantity is the 2-dimensional 
vector sum of the transverse momenta in C, 



,c 



(5-1) 



(with Q±,c defined in eq. (|4.1|) ) which, by conservation of momentum, is equal to 
(minus) the vector sum of the transverse momenta outside C. Thus we obtain a 
whole series of 'indirectly global' event shapes, whose definitions are similar to those 
of section |L but with £g replaced by TZ± t c- 
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(5.2a) 
(5.2b) 
(5.2c) 
(5.2d) 
(5.2e) 



Experimentally, we envisage that the main difficulty that will arise specifically 
for this class of observables is the accurate determination of H± ! c, since it involves 
a cancellation between hard momenta. The extent to which this 'missing transverse 
momentum' can be well measured will determine the extent to which it will be 
possible to study the region of low event-shape values. 

The single-emission properties for two of these observables are illustrated in 
tables ^ and ||. As expected, the outgoing legs have identical properties to the global 
variants of the observables. The incoming legs have = a, and a simple calculation 



- 22 - 



leg t 




bt 


9i{4>) 


di 


1 


1 





1 


1/sine* 


2 


1 





1 


1/sine* 


3 


1 


1 


sin 2 (p 


1/sin 2 9* 


4 


1 


1 


sin 2 (p 
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Table 7: Leg properties for tj_ 
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1 
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sin0|)/2* 


2/ sin^* 
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1 
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sin0|)/2* 


2/ sin^* 
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1 







sin 


2/ sine* 


4 


1 







sin (j)\ 


2/ sine* 



Table 8: Leg properties for T m ^\ 
(ln 5l)2 (<^))=4G/7r-21n2. 



reveals that dg — 1/ sine and <%(</>) = 1 for most observables, including tl^. An 
exception is T m ^, for which there is an interplay between recoil dependence present 
(implicitly) in T m C , and the additional TZ^c term. 
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Figure 11: for r^. 
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Figure 12: for T m ^. 



The J-'(R') functions for t±k and T m ^ are shown in figures |TT| and [TJ. For 
small values of R', J-'(R') is below 1, as for the observables discussed in the previ- 
ous sections. However for sufficiently large R', F(R') starts to increase and it has 
a divergence at some finite value of R' = R' c . Such divergences are actually a well- 
understood phenomenon [^, ^] , characteristic of observables for which contributions 
from multiple emissions can cancel. Schematically, they occur because for sufficiently 
large R' it becomes more favourable to suppress the value of the observable via can- 
cellations between emissions than by Sudakov suppression, and the resulting change 
in parametric behaviour of E-^(t») cannot be represented as a subleading correction 
to a Sudakov resummation. 

For the observables being discussed here, the cancellation is in the vector sum 
of transverse momenta. This being a two-dimensional cancellation, the divergence 
is located at R[ +2 = 2 [19|, |65|], where R[ +2 is the part of R! associated with the 
incoming legs. For cases where R[ +2 = hR', such as the gg — ► gg channel for T m ^ 
then R' = 4. For the gg — > qq scattering channel, in which R' comes mostly from 
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the incoming legs, the divergence is at lower R'. Furthermore for any given fixed 
channel, the divergence occurs earlier for the transverse thrust than for the minor, 
because the thrust has smaller LL terms, and accordingly a larger proportion of R' 
is associated with the incoming legs. 

The divergence is a structure that appears even more strongly at yet higher 
orders (NNLL, . . . ) and it could in principle be resummed. Technically speaking, one 
would have to carry out a 6-space resummation |68[] on the incoming legs, a 'normal' 



resummation for the outgoing legs, and include a non-global resummation j37| to deal 



with the dynamic discontinuous non-globalness [^0[ associated with the boundary of 
C. This is technically rather challenging and so far the only such resummation [p5|1 , 
for a jet-broadening in DIS, neglected the non-global logarithms. The conclusion 
of that study (for a case in which R' c = 4) was that the further resummation had 
negligible practical impact except in a region where the distribution was already very 
strongly suppressed. 
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Figure 13: The distribution of t± h for 
separate hard-scattering channels. 
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Figure 14: The distribution of T m ^i for 
separate hard-scattering channels. 



We therefore follow the approach of |65j and consider the distributions for the 
observables without any additional resummation, figs. [13] and |14[ The distributions 
are cut (somewhat arbitrarily) at the point where R' = 7/8R' c — one sees that 
this corresponds to different values of lnt> , depending on the resummation channel. 
However for each channel most of the cross section is in the theoretically reliable 
region where the divergence can be ignored. We note that some problems do arise in 
the combination of channels with divergences at different points, notably from the 
point of view of establishing the region in which the sum over all channels is well 
estimated. The detailed discussion of this issue will be left to future work f53| . 

We conclude this section by observing that there are reasons to believe that 
these observables, though global, should have a relatively limited sensitivity to the 
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underlying event. This is because in the sum of perturbative and 'underlying-event' 
contributions to 7l±,c, after integration over azimuthal angles, the leading piece of the 
underlying event component averages to zero, leaving only a subleading component. 12 



6. Conclusions 

In this article we have defined and resummed a set of dijet event-shape observables 
and jet rates that are suitable for study at hadronic colliders, both from a theoretical 
and an experimental point of view. The principle issues that had to be reconciled 
were those of globalness, which greatly simplifies the resummed theoretical predic- 
tions, and of limited experimental reach in rapidity. At first sight these requirements 
seem contradictory, however we saw that it is possible to define several classes of 
observables which meet both criteria. 

For a majority of the observables, the distribution is concentrated in a region 
where the limited experimental rapidity reach can be ignored in the resummed pre- 
diction. This is the case for nearly all the directly global observables of section |3] 
and those with exponentially suppressed forward terms, section f|, the latter having 
been specifically designed to be optimal in this respect. The observables of section D, 
defined only in terms of particles in the central region, avoid the problem altogether, 
implementing globalness 'indirectly' via a recoil term (though this affects the validity 
of the resummation in a limited region of very small event-shape values). 

The study of a range of observables, as presented here, would have been consid- 
erably more difficult without the help of the automated resummation tool, CAESAR. 



For brevity we omitted a range of further results [62]. Full phenomenological studies 

some 



will necessitate also matching with fixed order calculations such as p3| , \L% p~5 
relevant issues for 'practitioners' having been discussed in appendix [B]. 

One of the main reasons for examining several classes of observable (and a num- 
ber of observables within each class) is that they have quite complementary sensi- 
tivities to different physics issues, both perturbative and non-perturbative. Possible 
studies include measurements of a$, tests of novel perturbative QCD colour evolution 
structures that arise in events with 4-jet topology, and investigations into hadronisa- 
tion effects and the properties of the underlying event. Overall, we believe that this 
represents a potentially far richer programme of investigations than has, for example, 
been carried out in e + e~ or DIS processes. 

Acknowledgements This paper has evolved from discussions with a number of 
colleagues, including in particular Leonard Christofek, Yuri Dokshitzer, Joey Huston, 



12 Similar arguments have been applied for non-perturbative studies of jet-broadenings in e + e 



and DIS |65|, p9| . with a good degree of phenomenological success |40| [70|, |7lj, [72[. The argument 
assumes a lack of azimuthal correlation between the underlying event and perturbative radiation 
from the 'hard' event. 



- 25 - 



Pino Marchesini, Aurore Savoy-Navarro and Mike Seymour. We are grateful also the 
CERN and the IPPP, Durham for the use computing facilities. GPS would like to 
thank also the Fermilab theory group for hospitality while this article was being 
prepared. 



A. Scale choices 



As discussed in [[H], |65|, [73| , in a resummation at NLL accuracy there is an intrinsic 
ambiguity in the choice of the logarithm to be resummed, for instance instead of 
resumming 

L = ln-, (A.l) 

v 

one could equally well want to resum 

L = ln^-, with X = 0{1). (A.2) 

This alternative choice for the logarithm affects the functional form of the single- 
logarithmic terms, but not the overall answer at NLL. Changing X allows one then 
to estimate the size of a specific class of higher order (NNLL) corrections. In two-jet 
event shapes, the default value for X was chosen in a 'natural' way by requiring that 
the coefficient Gu, the O (a s ) single logarithmic term, contained only terms due to 
hard collinear splitting. This corresponded to taking X = 1/d (usually d\ = d 2 = d 
in e + e~) and ensured that one didn't introduce large spurious subleading logarithms 
whose coefficients involved powers of In d, associated just with the overall (arbitrary) 
normalisation of the observable. 

The equivalent procedure more generally would be to define X so as to cancel 
the coefficient dege((p) in eq. (|2.6|). The obvious difficulty though is that dege(<fi) is 
not a constant — it depends on the hard leg i and on the value of the azimuthal 
angle. 



Inspecting the detailed resummation formulae, in particular eq. (3.6) of [BJ, one 
sees that the problem of the azimuthal dependence can be eliminated, insofar as the 
resummed result depends only on In dg = In di + J # ln^(0). Furthermore, for each 
leg (if all the bi values are the same), the term involving lno^ comes in with a weight 
proportional to the colour factor, C#, of the leg. So for X to cancel the coefficient 
dege(<f>) on average, it suffices to take 

-. n n 

\nX = -—Y j (C l \ndt) , C T = y^C e . (A.3) 

Ct 1=1 e=i 

with n the number of hard legs. This choice means however that L is defined dif- 
ferently for different colour configurations, which while legitimate, may seem a little 
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unnatural. In this study we therefore adopt a prescription for X, which is numer- 
ically almost indistinguishable from the choice in eq. ( |A.3|) , but does not have the 



unusual feature that the logarithm depends on the specific colour configuration, 

1 - 

lnX = Vlnd £ . (A.4) 

n 

1=1 

For situations in which all legs have the same colour factor, it is equivalent to 
eq. ([Op. 



B. Comparison to fixed order 

A valuable cross-check on any resummed prediction is that its expansion to fixed 
order coincides with the logarithmically enhanced terms of exact fixed-order predic- 
tions. For e + e~ — > 2 jets and DIS 1 + 1 jet resummations, such checks have been 
carried out to NLO accuracy. The current status of the NLO codes suitable for 



NLO prediction of hadronic dijet event shapes is that NLOJET++ is available 
publicly, while TRIRAD [|15j can be requested from the authors. 

We have carried out comparisons with NLOJET++. 13 This turned out to be 
feasible at LO (for the a s L 2 and a s L terms of S^(v)/cj^), while we encountered dif- 
ficulties in obtaining a determination of the NLO contribution that was sufficiently 
accurate to enable a meaningful comparison with the expansion of the resummations. 
Nevertheless even the LO comparisons represent a non-trivial check of the resumma- 
tions, in particular in view of the involved four-jet large-angle colour structure. 

Fixed-order calculations are relevant not only as checks of the resummation, but 



more importantly for matching [|3£|, so as to obtain predictions that are valid for 
large as well as small values of the observable. The matching also gives a partial 
improvement in the accuracy at small values of the observable, supplementing, for 
instance, eq. ( |2.5| ) with a term known as C±, which depends on the momentum 
configuration and the channel, 

fB,s{ v ) = (l + Cl > B > 5 ~^ + • • •) ex P [Lgi,s{oisL) + 92,BA a sL) + a s g 3! B,s(a s L) H ] , 

(B.l) 

where we have explicitly indicated the dependence of the g n (a s L) on the channel 5 
and (for n > 2) on the momentum configuration B (elsewhere this dependence has 
been left implicit). 

The improvement in accuracy comes, for example, because the product of Ci,b )( $ 
and the first term in the expansion of Lgi^(a s L) helps to fix the a™L 2n ~ 2 terms in the 
expansion of fjs t s{v), or equivalently, after summing over channels and integrating 



13 We are grateful to Zoltan Nagy for having provided us with a prerelease of a new version of his 
code. 
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over Born configurations, in E^(v)/cr^. Currently however, fixed-order programs are 
at best able to provide a weighted sum of Ci^s values across all channels 



(Ci, B )i 



6 dB 



(B.2) 



Since different channels have different colour factors appearing in Lgij(a s L), the 
averaging over channels for C\^,& means that one cannot reconstruct the full infor- 
mation for the sum of products of C\^,s and Lgi t s(ce s L) 



C liB Lg 1 (a s L)) 5 ^ (C 1;B )s {Lgi{a s L)) s . 



(B-3) 



The fact that the Ci B s are characteristics of the soft and collinear limit, means that 



they can in principle be unambiguously extracted channel-by-channel |xj , enabling a 
proper average to be carried out. However the identification of the channel requires 
that one have information on the flavour of individual partons in the fixed-order 
calculation. Though such information is present in some form in existing fixed-order 
codes [|3|, pL5| , it is not, as far as we understand, available through the external 'user' 



interfaces to those codes. The availability of a method to access this information 
more straightforwardly would be of considerable help for the matching. 
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